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The viscosity of liquids under nanoconfinement remains controversial. Reports range from spon- 
taneous solidification to no change in the viscosity at all. Here, we present thorough measurements 
with a small-amplitude linear atomic force microscopy technique and careful consideration of the 
confinement geometry, to show that in a weakly interacting liquid, average viscosity remains bulk 
like, except for strong shear thinning once the liquid is confined to less than four molecular lay- 
ers. Overlaid over this bulk-like viscous behavior are stiffness and damping oscillations, indicating 
non-continuum behavior, as well as an elastic response when the liquid is allowed to order in the 
confinement gap. 

PACS numbers: 68.08.-p, 66.20.-d, 07.79. Lh, 62.10.-(-s, 61.30.Hn 



Nanoconfined liquids play an important role in 
nanotechnology[l[ tribology [3|], geology [J], biology 
and medical technology [g, but their properties re- 
main controversial Although it has been widely 
accepted that nanoconfined liquids are different from 
the bulk, their often surprising viscoelastic behavior is 
not yet completely understood. According to some ex- 
perimental groups, liquids pass through a spontaneous 
phase transformation from liquid to solid below a certain 
threshold confinement • Other groups see an enhance- 
ment in viscosity or relaxation time only [9|_jOr see very 
little change in viscosity upon confinement 19[. The be- 
havior may also be different for different liquids and de- 
pend on the dynamic rate of confinement IB" IM IS 21 1 ■ 

Atomic force microscopy (AFM) or surface force appa- 
ratus (SFA) measurements of the mechanics of nanocon- 
fined liquids either measure a damping coefficient under 
compression, involving squeeze-out [2 2], or employ shear 
measurements [1, [13] • While the latter leads to a more 
straight-forward measure of effective viscosity, squeeze 
out measurements are often easier to perform and can 
be combined with measurements of the compressive elas- 
ticity of the nanoconfined liquid. Moreover, some recent 
reports suggest that the behavior of the liquid may be 
anisotropic, i.e. viscous behavior may be different de- 
pendent on the load direction (2l| . 

A problem with comparing measurements is assessing 
the influence of confinement size. In order to normal- 
ize the compression force, the Derjaguin approximation 
is commonly usedi23|. For a sphere-on-fiat, this approx- 
imation states that the quantity F/R is independent of 
confinement geometry, where F is the force and R is the 
radius of the sphere, which in AFM is the cantilever tip. 
The Derjaguin approximation holds over many orders of 
magnitude, and allows for the comparison of SFA and 
AFM measurements in a variety of geometries. For a 
sphere-on-flat, the interfacial energy per unit area can 
be directly related to the interfacial force, as we have 



w = dW/dA = -l/{2TTR)dW/dh = F/2TrR, where h is 
the vertical distance along the sphere, measured from the 
substrate to the sphere surface. Since the interfacial en- 
ergy is assumed to be independent of radius, F/R should 
be a the same for any radius sphere approaching a flat 
substrate. We use a small amplitude dynamic technique, 
which does not measure force, but rather the local force 
gradient or interaction stiffness (where kmt = —dF/dz) 
and the squeeze damping, C, allowing determination of 
the effective viscosity. For the interaction stiffness, we 
find k = -dF/dz = d^W/dh^ = -d/dh{2TrRdW/dA) = 
—2TrRdw/dh. Assuming that the gradient of the interfa- 
cial energy is a function of h only, the interaction stiffness 
k should also scale linearly with the radius of the sphere. 
Note that typically dw/dh < 0, and k is positive. 

The correct expression for the squeeze damping term 
can be found by classical lubrication theory from 
Reynolds equation, which, in the axisymmetric case 
reads [H: 



19/ rh^ dp 
r dr \ rj dr 



-12 



dh 
'~dt 



(1) 



For a sphere with h{r) = Hq + R — \/R^ — r-^, constant 
dh/dt — V and constant r/, ([T|) can be integrated twice to 
give 



F = Cv^ 
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Typically, only the leading term is used, as the correction 
terms become very small for Hq <C R- However, equation 
([T]) can also be integrated for an arbitrary axisymmetric 
shape of the tip (for non-symmetric tips, the more general 
two-dimensional form of the Reynolds equation must be 
solved, which is significantly more involved) and for non- 
constant viscosities. 

Measurements were performed on nanoconfined 
tetrakis(2-ethylhexoxy)silane (TEHOS), an inert liquid. 
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FIG. 1. Scanning Electron images of the four probes used in 
this study. 



which has been well characterized through X-ray [25j, flu- 
orescence correlation [26| and AFM measurements 13 1 . 
Stiffness and damping were determined with a sub- 
resonance dynamic AFM technique, by measuring am- 
plitude and phase of the approaching tips at amplitudes 
of < 0.6 nm and an approach speed of 1 nm/s. The AFM 
cantilevers were completely immersed in the liquid dur- 
ing the measurement avoiding any surface tension effects. 
A home-built, fiber-interferometric AFM was used for 
the measurements 27|. Cantilevers, substrates and the 
measured liquid were thoroughly cleaned prior to mea- 
surements. We used four different cantilevers, with stiff- 
nesses 1.0,2.4,35 and 55N/m, with nominal tip sizes of 
0.5,1.4,4.7 and 18/^m, respectively, in order to explore 
the effect of tip size on the stiffness and damping of 
nanoconfined TEHOS. 

Figure [T] shows scanning electron microscopy (SEM) 
images of the four tips used in this study. Such imaging 
is limited, but allows us to approximate the tip shape 
sufficiently to get close matches between the measured 
and calculated damping curves. As seen, the radius of 
curvature Rc of the various tips is not necessarily the 
same as the nominal radius of the tip end, R. This 
is due to flattening of the tip during the initial mea- 
surement stages. This flattening has the advantage that 
small asperities are eliminated, and unlike other reported 
measurements [28l |. our measurements reproduce the ex- 
pected Derjaguin dependence of the stiffness on the ra- 
dius quite well (Figure [2]) . We found that the stiffness 
is linearly dependent on R and not on the radius of cur- 
vature of the tip, Rc- This is not too surprising, as the 
Derjaguin approximation relates to the contact radius, 
which is mostly determined by actual tip size. 

Figure [3] shows stiffness and damping for the largest 
probe. Stiffness and damping reflect molecular layering 
of the liquid, exhibiting regular peaks roughly commen- 
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FIG. 2. Reduced stiffness k/R for the largest R = IS/^m and 
the smallest tip radius R = 0.5/^m. Arrows indicate layering 
of the liquid, with a temporary pinning of the tip and local 
peaks in the stiffness. Unlike our previous reports, we were 
able to extract the actual tip position and therefore could plot 
these data against the actual tip-surface distance. Dashed 
arrow shows an intermediate peak apparent in both data sets. 
Inset: Log- log plot of stiffness k versus tip radius R for all 
four tips. Dashed line has slope of one, indicating a linear 
relationship between k and R. Triangles: h = ?, nm. Circles: 
h = 2 nm. Squares: h = 1 nm. 



surate with molecular size. It should be noted that seeing 
these peaks with a colloidal probe tip is challenging, as 
any tip roughness will destroy average ordering in the 
confined liquid. To our knowledge, this is the first time 
that stiffness and damping oscillations in a nanoconfined 
liquid have been seen using a large colloidal AFM probe. 
The stiffness peaks are at locations where the tip remains 
stationary, i.e. where the tip is pushing against an inte- 
ger number of molecular layers. By contrast, the damp- 
ing peaks preferentially line up between the stationary 
points (dashed arrows), indicating that enhanced damp- 
ing may result from rearranging molecules as the liquid 
transitions from iV to — 1 layers between tip and sub- 
strate. However, as shown previously, the fact that stiff- 
ness and damping peaks are out of phase with each other 
is not universal, but depends on the squeeze rate 17. isj. 
In TEHOS, we almost exclusively observe the damping 
and stiffness to be out of phase, indicating an elastic-like 
behavior of the ordered molecular layers of the liquid. 

Figure U] shows the damping curves measured with the 
four probes on a log-log plot. Generally, the damping is 
proportional to h^^ as expected, although there are de- 
viations. At small separations of a few nanometers, and 
therefore a few molecular layers, damping saturates. This 
saturation of the damping can have two origins. Elasto- 
hydrodynamic effects can lead to a deformation of the 
substrate by the localized liquid pressure above it[29|. 
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FIG. 3. Stiffness and damping for the largest {R = 18^m) 
probe. Dotted lines indicate layering of the liquid. Arrows 
point out increased damping between the layers. 



However, if we calculate the critical liquid gap distance 
where such effects should be seen, it is of order 10~^^ m, 
which is much too low to account for the saturation seen 
at a few nanometers in these measurements. 

Alternatively, a saturation in the damping may sug- 
gest a reduction of the effective viscosity at close ap- 
proach. Since at close approach the squeeze-out leads to 
increased shear in the liquid, this suggests shear thinning. 
This is consistent with previous reports on shear thin- 
ning effects in nano confined octa-methyl cyclotetrasilox- 
ane (OMCTS) ^ and in water [l6|^hear thinning can 



be described by the Carreau model[3l|, where the effec- 
tive viscosity is given by 



V = 



1 + (A7) 



l-s 



(3) 



where A is the characteristic time scale at which shear 
thinning becomes observable, and s is the so-called shear 
thinning exponent (for shear thinning, s < l).The shear 
rate can be estimated by calculating the volumetric flux 
per unit area as the tip squeezes out the liquid. We find 
j = Ldy. — ^L-^ = ™. The average shear rate at an 
axial distance r is then the flux divided by film thickness, 
or 7 = This can be substituted into equations (P) 
and Qwhich can then be numerically integrated for any 
tip shape h = h(r). 

We estimated the tip shape from the SEM images, from 
which a reasonable range of possible radii of curvature 
can be obtained. Within this range, we adjusted param- 
eters to obtain the best fit to the measured data. In this 
sense the tip shape is consistently determined by both 
imaging and fitting to the obtained damping curve. It 
should be noted that at large distances, before the damp- 
ing saturates, s and A have practically no influence on the 
fit. Thus the fit, together with the SEM imaging, deter- 
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FIG. 4. Log-log plots of the damping coefficient of the con- 
fined liquid for all four tips. Tip parameters are given in the 
respective figures. Dashed line indicates slope without assum- 
ing any shear thinning. This slope is typically equal to -1, as 
expected for a spherical tip, but deviations can be observed 
if significant flattening is present. 



mines the effective tip geometry quite well. We also kept 
the fitting parameters A and s as close as possible across 
all measurements. We found that a value of s = —0.7 
gives best fits for all tip geometries used, while A varied 
in the relatively narrow range of 90-180 ms. For all fits, 
we used a viscosity of 0.009 Pa s, which is the expected 
bulk viscosity of TEHOS at 293 K [s^. This univer- 
sality of our fits to a variety of tip geometries gives us 
confidence that no viscosity enhancement was observed 
in these measurements. The saturation of damping at 
close approach suggests shear thinning, which is consis- 
tent with previous reports 



16,aQl. It should be noted 



that the fact that shear thinning is observed when the 
closest tip-surface distance is just a few nanometers does 
not mean that shear thinning occurs at this separation. 
The tip is curved, and therefore there is a range of liquid 
thicknesses at different locations along the radius of the 
tip. Moreover, maximum shear fiow does not occur at 
the tip center. 

In conclusion, we have shown that in a simple, weakly 
interacting liquid, the average viscosity measured by 
squeeze-out is bulk-like even for molecularly thin liquid 
layers. However, at moderate shear rates and very close 
to the surface, shear thinning is observed. Overlaid on 
the average bulk-like/shear thinning behavior are oscil- 
lations in viscosity commensurate with molecular layer- 
ing of the liquid. The fact that damping and stiffness 
peaks are out-of-phase with respect to molecular layer- 
ing indicates dynamic effects, which lead to an elastic 
response when the liquid is molecular ordered between 
the two confining surfaces 17|, |l8| . Thus the behavior of 



nanoconfined liquids can be surprisingly rich, exhibiting 
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average continuum behavior, down to molecular scale, 
which is however modulated by discreteness of the sys- 
tem, once it is confined to less than 4-5 molecular lay- 
ers. This is also reminiscent of plastically deformed solid 
systems, which can show combine discrete jumps upon 
deformation, superimposed on an overall continuum-like 
response j33{. 
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